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We examine a parametric cycle in the A^-body Lieb-Liniger model that starts from the free system 
and goes through Tonks- Girardeau and super- Tonks- Girardeau regimes and comes back to the free 
system. We show the existence of exotic quantum holonomy, whose detailed workings are analysed 
with the specific sample of two- and three body systems. The classification of eigenstates based on 
clustering structure naturally emerges from the analysis. 
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I. INTRODUCTION 

Among the solvable models of quantum mechan- 
ics, Lieb-Liniger system [l| belongs to the selective 
class of models that are genuinely many-body. It 
is a system made up of identical bosons interacting 
through two-body contact force. It has later been 
shown that the one-dimensional system of identical 
fermions with two-body contact interactions can be 
rigorously mapped to Lieb-Liniger system with strong 
and weak coupling regimes interchanged [2, y] . Several 
further extensions of the model with anyon statistics 
has been found [^-|2|; and they are also known to be 
mathematically equivalent to original model. Ther- 
modynamics of Lieb-Liniger model has been exten- 
sively studied |8l4ll|. 

What has made the Lieb-Liniger model a focus 
of renewed recent attention is its experimental re- 
alisation in the form of Tonks- Girardeau gas 12- 
ll4J . It has been shown that the coupling strength of 
the Lieb-Liniger system can be experimentally con- 
trolled through the Feshbach resonance mechanism 



15|. In recent experiments by Haller and collabora- 



tors [16|, |l7[ : a smooth change of the coupling strength 
from large negative values to large positive values, 
where one finds the super Tonks-Girardeau system 
18|, has been realised. 
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The continuous transition from a strongly repul- 
sive to strongly attractive regimes of Lieb-Liniger 
model inspires us to propose following parametric cy- 
cle C: We start with the non-interacting limit, in- 
crease the coupling strength adiabatically, reaches 
strongly attractive regime crossing the ±oo coupling 
limit, then decreases the absolute value of negative 
coupling strength until it reaches the non-interacting 
limit again. In this paper, we show that the initial en- 
ergy eigenstates of the cycle are different from the final 
eigenstates, although the initial and the final Hamil- 
tonians are identical. 

This phenomenon, the so-called exotic quantum 
holonomy, in which quantum eigenvalues and eigen- 
states do not come back to the original ones af- 
ter a cyclic parameter variation 19[, belongs to a 
wider class of quantum holonomy that comprises both 
the celebrated Berry phase [20| and the Wilczek-Zee 
holonomy 21| which appears in systems with degen- 
erate eigenvalues. The exotic quantum holonomy in 
(5-function potential system has been already consid- 
ered in 123. Here we report a novel finding of the 
quantum holonomy in many-body systems interacting 
through (5-function potential. 

The plan of this paper is as follows: In Sec. |lll we 
derive the spectral equation for Lieb-Liniger model 
in two different forms to demonstrate the presence of 
quantum holonomies with respect to C. In Sec. lIIIi we 
show that the backward cycle is not always possible 
due to the clustering of particles. This leads to the 
concept of minimal states, which we utilize to classify 
the spectrum of the system in Sec. IIVI In Sec. |Vl we 
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provide another view of the quantum anholonomy by 
focusing on the two-body system through the com- 
plexification the couphng strength. Sec. |Vl] contains 
our conchision. 



II. ADIABATIC CYCLE C FOR 
LIEB-LINIGER MODEL 

Let us consider N bosons confined in a one- 
dimensional space. The system is described by the 
Hamiltonian 



1^ 92 



N j-1 



where the unit is chosen such that h and the mass of 
a particle can be set to unity. The parameter g is the 
interaction strength. We impose the periodic bound- 
ary condition to the position space with the period 
27r. 

We look at the dependence of eigenenergies and 
eigenvectors on the coupling strength g. In partic- 
ular, we focus on the cycle C, which consists of three 
stages C^"^ (s = 1,2,3): In the first stage C^^'> , g is 
prepared to be and, is adiabatically increased to 00; 
Next, in the stage C^'^\ g is suddenly flipped from cxd 
to —00; In the final stage C^^\ g is again adiabatically 
increased to 0, which is the initial value of g. We de- 
note the initial and final points of C as g = and 
g = 0—, respectively, to distinguish them. 

The eigenvalue problem of H, Eq. ([T]) can be solved 
by the Bethe ansatz, where an eigenfunction is com- 
posed by N plane waves specified by a set of quasi- 
momenta, also called rapidity kj, that satisfy 

kji + \g 



exp (27rifcj ) = 1 [ 
1^3 



kji ~ iff ' 



(2) 



a. 



where kij = ki -- n,^ 

We examine how kj{gys, which are chosen to be 
smooth as g is varied, are changed by the cycle C. 
The function kj (g) completely characterizes the para- 
metric evolution of eigenenergies as well as the "adia- 
batic" evolution of eigenvectors along C. The analysis 
is decomposed into the three stages C^'*) (s = 1, 2, 3). 

At the initial point 5 = of the first stage C^^\ kj{0) 
takes an integer value. Without loss of generality, we 



can choose the order of kj (g) 's so as to satisfy fci {g) < 
^2(5) < ••• < k^ig) for small positive g [23|. This 
ensures fci(O) < fc2(0) < • • • < kN{Q) 

We introduce two quantized quantities which is con- 
served during the parametric evolution of kj {g) along 
C^^\ Such "topological invariants" provide a way to 
evaluate the change of kj (g) induced by the stage C^^^ . 

Firstly, during the interval < g < 00, we have an 
integer 

1 



Ijig) = kj{g) 



Eg 
arctan -— . 



(3) 



This is a consequence of Eq. ([2]) and (i + i)/(t — i) = 
_g-2iarctant^ which is applicable as long as f^^ 7^ 0. 
We use the principal branch of arctan throughout 
this manuscript. This is justified for Eq. ^ because 
g/kji{g) does not across its standard branch cuts. 



which emanate from ±?' to ±ioo [24|. The continu- 
ity and discreteness of Ij (g) in < .9 < 00 imply that 
Ij{g) takes a constant value, which can be determined 
from the limit p J, 0, i.e.. 



I,{g) = k,{Q). 



(4) 



Secondly, Eq. ([2|) and another formula for arctan 
implies that, in the interval < 9 < 00, 

Jj (5) = kj (5) + - X! arctan — 



(5) 



is a half-integer for even N and an integer for odd 
N [l|. Following a similar argument for Ij{g) above, 
we obtain the value of the invariant Jj {g) for < g < 



•^■(5) = kjioo). 



(6) 



Now we evaluate the change of kj{g) during C^^^ 
using these invariants. From Eqs. ([3]) and (O, we 
obtain 

fc,(oo)-fc,(0) = i^sgn5R^^, 



(7) 



lii3 



where we used the identity 



arctan(f) + arctan(l/t) = — sgn[3ff(t)]. (8) 

We note that the right hand side of Eq. ([7]) makes 
sense only for < .g < 00. Here, kji(g) is positive for 
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j > I and negative for j < I, since we have assumed 
the order of kj{g) at the initial point of C^^\ and the 
sign of kji{g) docs not change for g > [23|. This 
imphes 

.kjiig) 



9 



yj sgn3?- 
Accordingly, we obtain 



1=3 + 1 



1 = 1 



kjioo) -kj{0) = j - 



N + l 



(9) 



(10) 



Next we examine the second stage C'^^\ where g sud- 
denly changes from oo to — cxo. Note that all A;j(oo)'s 
are finite because of Eq. ([TO)). Since a finite root of 
the Bethc equation, Eq. ([2]), at g = cxd is also its root 
at g ~ —CXD, we employ a smooth extension of kj{g) 
along C^^\ i.e., 

fc,(-oo) = fc,(oo). (11) 

The analysis for the final stage C'"^-' resembles to the 
one for C^^\ First, Jj{g) is a half-integer for even TV 
and an integer for odd A^ in the interval — cxd < g < 0. 
In order to ensure that Jj{g) is an invariant within 
this interval, we need to assume that kji/g does not 
cross the branch cuts of arctan. The proof of this as- 
sumption is to be published elsewhere |25|. Here let 
us just note that it is true for N ~ 2, and is numeri- 
cally confirmed for a numbers of lower-lying states of 
N = 3 and 4 in the following sections. Second, Ij{g) is 
also an integer for — oo < .g < 0. We also assume that 
g/kji does not cross the branch cut of arctan in this 
interval. This implies kj{0—) = Ij{g)- The change of 
kj{g) in the path C^"^-* is given by 

k,iO-)-k,i-^)^-lj2sgn^'^^'^^^ 



1^3 



9 



(12) 



We can ensure that 



^Mg)] K^Mg)] <---<^[kN{g)], (13) 

because it holds at g = — oo. Recalling the fact that 
g is negative here, we obtain 

k,iO-)-kA-oo)^j-^. (14) 

Combining above three arguments, we obtain a non- 
trivial change of kj{g) due to C in the form 



Note that the total momentum remains unchanged 
during the cycle C. The final energy and state after 
the adiabatic cycle, however, are different from the 
initial ones, showing that C induces the eigenenergy 
and eigenspacc anholonomies [22| . We also remark 
that ki < k2 < ■ ■ ■ < kjsi holds at the end of C. 
This implies that we can repeat the adiabatic cycle 
C arbitrarily, and, the repetition of C will induce the 
further instances of the eigenenergy and eigcnspace 
anholonomies. 

We can summarize our results in terms of a mapping 
between two sets of quasi- momenta of free bosons, i.e., 
fcj(0)'s and fcj(0— )'s . It is sufficient to consider the 



case that initial condition riw 



fcj(O) satisfies rii 



< 



^2 < •■• < riN- With the notation n' = fcj(0— ), 



the mapping (ni, ^2, . . . , njv) ^-^ (n'l, 
F{ni, 712, ■ • ■ , 'T-Tv), which is given by 



'-2J- 



>"W) 



fc,(0-)-fc,(0)=2j-(iV + l). 



(15) 



F{ni,n2,...,nN) 

= (711 - iV + 1, 772 - A^ + 3, . . . , 7lAr + A - 1), (16) 

expresses the quantum holonomy induced by the cycle 
C. 

III. INVERSE CYCLE 

We now examine the inverse of the cycle C. In con- 
trast to the forward cycle C, the parametric variation 
along the inverse C~^ is not always possible. This is 
because the clustering of particles aX g — — oo induces 
the divergence of eigenenergy [1|. Such a clustering 
invalidates the use of the Hamiltonian, Eq. ([T]). We 
call an eigenstate of free boson at g = a minimal 
state if the the parametric variation along C^^ is im- 
possible. The precise condition for appearance of the 
minimal state is the subject of this section. 

Formally, C~^ corresponds to the inverse of the 
mapping F^ Eq. (|16p . on the sets of quasi- momenta 
at g = 0: 

i^"^(77l,772,...,777v) 
= (771 + Af - 1, 772 + A^ - 3, . . . , Tlw - A^ + 1), (17) 

where we impose the ordering condition 77i < 772 < 
• • ■ < 77 jv. When the distance between 77 jS arc far 
enough, F^^ preserves the ordering. This is the case 
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that C~^ can be realized, and the resultant energy 
and quantum state are the solution of the eigenvalue 
problem of H, Eq. ([T]), at g = 0. On the other hand, 
when a pair of rijS is too close, F^^ breaks the or- 
dering, which implies the emergence of the cluster- 
ing of particles during the inverse cycle. There are 
two possible cases. The first case is where a pair of 
quasi- momenta, say, Uj and rij+i are degenerate, i.e., 
Uj = Tij^i. By applying F^^, the resultant quasi- 
momcnta satisfy rij > rij^i. In fact, the eigenenergy 
diverges — cx) as g — > — oo during C~^. The second 
case, Uj — nj+i -f- 1, also leads the clustering of parti- 
cles. 

The argument above is sufhcient to determine the 
condition for the minimal states. When there is, at 
least, a pair of two quasi-momenta at g = that sat- 
isfies 



Wj -nj+i\ < 1, 



(18) 



states specified by Uj and rij-i-i arc minimal states. 



IV. CLASSIFICATION OF SPECTRA 

Because of the existence of quantum holonomy, 
some states are reachable by the repetitions of para- 
metric cycles C and C^^ starting from one particu- 
lar eigenstate, while other states are not. This offers 
the classification of whole eigenstates into families of 
states connected by quantum holonomy. Such a family 
can be specified by a minimal state introduced above, 
because an arbitrary eigenstate with a finite energy 
can become minimal by a finite repetition of C^^. 

From one minimal state, we can find other minimal 
states using the symmetries of the Hamiltonian ([T]). 
Suppose that a minimal state is specified by quasi- 
momenta (ni, 712, ... , tt-n). The translational symme- 
try implies that (ni -I- 1, n2 -f 1, . . . , njv + 1) is also 
a minimal state, whose total momentum is larger 
by N than the original one. For an arbitrary in- 
teger £, (ni + £,n2 + £^ . . . , n^ + £) is also a min- 
imal state. The refiection symmetry implies that 
(— n^r, . . . , — ri,2, — ni) is also a minimal state, which 
may or may not be different from the original state. 

Hence, it is sufficient to find all minimal states 




FIG. 1. Parametric evolution of eigenenergies of two- 
body Lieb-Liniger model, where x and y axes indicate 
(4/7r) arctan g and VE, respectively. The thick (black) 
and thin (blue) lines correspond to the families specified 
by the minimal states (0,0) and (0,1), respectively. See 
Eqs. (|2ip and (|22p . Note that the eigenenergies are con- 
tinuous at (? = ±oo. 



whose total momenta satisfy the condition 

TV v^ N 

2 ^-^ -^ - 2 



(19) 



to enumerate all minimal states using the translational 
symmetry, offering a way to classify the spectra of 
the Lieb-Liniger model completely. We illustrate this 
classification for a few-body cases. 

We start the analysis of A^ = 2 case with two mini- 
mal states 



(0,0) and (0,1). 



(20) 



We obtain two families of eigenstates at 5 = from 
these two minimal states, by repeating C 



and 



(0,0) H^ (-1,1) ^ (-2,2) 



(0,1) ^ (-1,2)h^ (-2,3) k^ 



(21) 



(22) 



respectively. The eigenenergies of these families are 
depicted in Fig. [TJ By shifting the total momentum 
from the two minimal states, Ec^. (1201) . we obtain an 
infinite number of minimal states {£,£) and {£,£ + 1) 
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FIG. 2. (Color online) Eigenenergies of TV = 3 families, 
where x and y axes are the same as in Fig. [I] The total 
momentum of all families shown here is zero. The thick 
(black) line corresponds to the family (Eq. (|24[) ) specified 
by the minimal state (0,0,0). The thin (blue) line cor- 
responds to the (—1,0, l)-family. These two families are 
trimers. The dotted (red), dashed (brown) and dashed- 
dotted (gray) lines are dimer families specified by mini- 
mal states (—1, —1, 2), (—2, —2, 4) and (—3, —3, 6) respec- 
tively. Although there are level crossings, the adiabatic 
theorem ensures that the adiabatic time evolution is con- 
fined within a family [331 ]. 



with an arbitrary integer £. The {i,i)- and {£,£ + 1)- 
families have the set of quasi-momenta at 9 = given 
by{(£-m,£ + m)}-^oand{(f-m,£+l + m)}-^o, 
respectively. This exhausts the minimal states and 
famihes for N ~ 2. 

The A^ = 3 case is far more complex than the N = 2 
case. First, we consider the case that the total mo- 
mentum is zero, where an infinite number of minimal 
states can be found. We depict some of them in Fig. [51 
There are two minimal states 



(0,0,0) and (-1,0,1), 



which are called trimers [26| , because the clustering of 
all three particles occurs in the limit g — > —00. The 
family of eigcnstates at g — specified by the minimal 
state (0,0,0) is: 



FIG. 3. (Color online) Parametric evolution of eigenener- 
gies of the N = 4: case, where x and y axes are the same as 
in Fig. [1] The thick (black) line corresponds to the family: 
(0,0,0,0) !-> (-3,-1,1,3) 1-^ (-6,-2,2,6).... The thin 
(blue) and dotted (red) lines correspond to (—1,0,0,1)- 
and (—1, —1, 1, l)-families, respectively. 



Besides, there are an infinite number of minimal states 
{(-^,-^,2£)},>o and {(-2£,^,^)},>o, (25) 

where the latter set can be induced through the use 
of the reflection symmetry. These minimal states are 
called dimers 



26[ , because the clustering of two parti- 
cles occurs in the limit g — > — cxd. Second, we consider 
the case ^jfij = 1. We have a trimer 

(0,0,1), (26) 

and an infinite number of dimers 

{(-£, -£, 2i + l)},>o, {(-^, ~l + 1, 2^)},>o, 
{(-2£+l,€,€)},>o, {(-2£,£,£-Kl)},>o. (27) 



Note that all minimal states that satisfy ^ 
v^"^J can be obtained from the minimal state with ^ 



-1 



j'h 



(0,0,0)k^(-2,0,2)k^(-4,0,4)^ 



(24) 



1 through the use of the reflection symmetry. We 
obtain all other minimal states from above using the 
translational and reflection symmetry. 

It is possible to enumerate minimal states and as- 
sociated spectral families in a similar way for larger 
A^. We simply close this section by showing several 
families of A^ = 4 system in Fig. [3l 
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V. EXCEPTIONAL POINTS 

So far we have focused on the quantum holonomy 
induced by the real cycle C. In Ref [27[, it is argued 
that, through an analysis of a quantum kicked top, the 
quantum holonomy has a correspondence with the ex- 
ceptional points, where the Hamiltonian exhibits non- 
Hermitian degeneracy 



28 



281. In other words, it is 



conjectured that the eigenenergy and eigenspace an- 
holonomy can be understood as a result of the meta- 
morphosis of eigenenergies and eigenstates induced by 
the encirclements around the exceptional points. In 
the following, we offer another example of this con- 
jecture using the two-body Lieb-Liniger model by de- 
forming C in the complexified (7-space. 

The analysis requires the non-Hermitian exten- 
sion of the Lieb-Liniger model [30|. We obtain 
eigenenergies with complex-valued coupling parame- 
ter g through numerical computation. We here fo- 
cus on (0,0)-family (Eq. dH]))- Let S„(.g) denote the 
eigenenergy of the state whose quasi-momenta take 
{—n,n) at g ~ Q. We depict En{g) for n = 0,1,2 
in Fig 131 We find that these eigenenergies compose 
a Riemann surface. Its Riemann sheets En{g) are 
connected by the exceptional points and associated 
branch cuts (cf. Ref. 



31 



Under the present choice of the branch cuts, all ex- 
ceptional points of (0,0)-family appear in the Eo{g) 
sheet. A pair of eigenenergies En{g) {n > 0) and 
Eolg) has a pair of degenerate points gn and 5*, where 
we choose 3(7„ < 0. We find that all degenerate points 
are of degree two. Hence the pair of eigenenergies for 
an exceptional point exhibits square-root type singu- 
larity. The encirclement around the exceptional point 
gn in complex g-plane induces the permutation of 
E(){g) and _E„(g). We numerically confirm these prop- 
erties of gn with n = 1, 2, . . . , 10. We find that 3?5„ 
and 3(7„ are decrease monotonically as n increases. 
We also obtain similar result for (0, l)-family. 

Let us consider the cycle that is a concatenation 
of C and Ci in Fig 2] (d) . Because this cycle encir- 
cles the exceptional point gi, the cyclic permutation 
{Eo,Ei) occurs. On the other hand, the cycle com- 
posed by C and C2 induces the cyclic permutation 
among {Eo,Ei,E2)- As the cycle involves more deeper 




-1 1 

Re(s) 



1 1 

Re(g) 



(c) 



(d) 



C 



% 




-t- 



-2-10 1 

Re(g) 

FIG. 4. (Color online) Contour plots of '^En{c): (a) n = 0; 
(b) 71 = 1; (c) 71 = 2. Lighter (darker) color indicates 
larger (smaller) value of 5R-E. Thin lines are the contours 
of 5R-E„(c). The exceptional points are indicated by filled 
circles. Bold lines indicate the branch cuts. While all 
complex EPs appear in -Bo(c), each En{c) (n = 1,2) has 
a single EP. (d) Schematic explanation of complex cycles 
that enclose exceptional points. We depict C by a thick 
(red) line. Dashed (blue) and dotted (black) curves indi- 
cate Ci and C2, respectively. See the main text. 

exceptional points, the accuracy of the the resultant 
permutation become better to approximate a shift to 
eigenenergies (£'0, -Ei, . . . ) n- {Ei,E2,. ■ ■), which is 
realized by the quantum holonomy along the cycle C. 
In this sense, we may say that the spectrum of Lieb- 
Liniger model feels the exceptional points that reside 
in the complex parameter space to induce the quan- 
tum holonomy along C. 

VI. CONCLUSION 

We have shown, in this work, that an eigenstate of 
the free Lieb-Liniger system g ~ +0 is transformed to 
another eigenstate with higher energy in the process of 
eigenspace anholonomy involving the parametric cycle 
g : +0 — > +CXD : —00 — > —0. Experimental tes ting 
should be within the range of current techniques [16 , 
llTJ . On the way to prove the existence of quantum 
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holonomy, we have demonstrated that the eigenstates 
of the Lieb-Liniger model can be classified according 
to their clustering property. The two and three boson 
system have been analysed in detail. 

Our result can be interpreted in terms of geometry. 
Consisting of real numbers and ±00, the parameter 
space of coupling strength is homeomorphic to S^. 
Therefore, our anholonomy is affected by the topol- 
ogy of S^. The presence of two kinds of invariants 
Ij{g), Eq. ([3|), and Jj{g), Eq. ([5|), for the parametric 
evolution of kj{g) reflects the fact that at least two 
charts are required for S^. Converting one of spec- 
trum condition to the other by using the formula of 
arctan ([5]) corresponds to coordinate transformation. 



The cycle of winding number m, C", increases fcj(O) 
hy m[2j ~ (N + I)]. 

The topological nature of the the quantum holon- 
omy implies that it is stable against, at least, small 
perturbations [321 ■ This also suggests that an experi- 
mental realization of the quantum anholonomy is pos- 
sible in one-dimensional bosonic systems. 
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